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Abstract 

We present measurements of the bb production cross section and angular 
correlations using the D0 detector at the Fermilab Tevatron pp Collider op- 
erating at y/s = 1.8 TeV. The b quark production cross section for \y b \ < 1.0 
and p\ > 6 GeV/c is extracted from single muon and dimuon data samples. 
The results agree in shape with the next-todeading order QCD calculation 
of heavy flavor production but are greater than the central values of these 
predictions. The angular correlations between b and b quarks, measured from 
the azimuthal opening angle between their decay muons, also agree in shape 
with the next-todeading order QCD prediction. 

PACS numbers: 13.65.Fy, 12.38.Qk, 13.85.Ni, 13.85.Qk - 14.65.Hq, 13.85.Qk 
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Measurements of the b quark production cross section and bb correlations in pp collisions 
provide an important test of perturbative quantum chromodynamics (QCD) at next-to- 
leading order (NLO). The measured b quark production cross section at \/s = 1.8 TeV [l|-|4j 
is systematically larger than the central values of the NLO QCD predictions |5|||. 

Measurements of bb correlations such as the azimuthal opening angle between b and b 
quarks allow additional details of b quark production to be tested since these quantities 
are sensitive to the relative contributions of different production mechanisms to the total 
cross section. Two measurements of bb angular correlations using the azimuthal opening 
angle of muons from the heavy quark decays, one at y/s = 1.8 TeV |7[] and another at \fs 
= 630 GeV H, are in qualitative agreement with the NLO QCD predictions. A different 
measurement at yfs = 1.8 TeV using the azimuthal opening angle between a muon from 
B meson decay and the b jet shows [|J qualitative differences with the predictions, while a 
direct measurement of bb rapidity correlations is found to be in agreement with the NLO 
QCD predictions. 

In this paper we provide an additional measurement of the b quark production cross 
section and bb angular correlations. The analysis makes use of the fact that the semileptonic 
decay of a b quark results in a lepton (here a muon) associated with a jet. We use a sample 
of dimuons and their associated jets to tag both b and b quarks. By tagging both the b 
and b quarks, we are able to significantly reduce the number of background events in our 
data sample. Also included is a revised measurement of the b quark production cross section 
based on an earlier DO analysis f|| using the inclusive single muon measurement. 



The DO detector and trigger system are described in detail elsewhere ||10|| . The central 
muon system consists of three layers of proportional drift tubes and a magnetized iron toroid 
located between the first two layers. The muon detectors provide a measurement of the muon 
momentum with a resolution parameterized by S(l/p) / (1/p) = 0.18(p — 2)/p©0.008p, with p 
in GeV/c. The calorimeter is used to measure both the minimum ionizing energy associated 
with the muon track and the electromagnetic and hadronic activity associated with heavy 
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quark decay. The total thickness of the calorimeter plus toroid in the central region varies 
from 13 to 15 interaction lengths, which reduces the hadronic punchthrough in the muon 
system to less than 0.5% of low transverse momentum muons from all sources. The energy 



resolution for jets is approximately 80%/ JE(GeV). 

The data used in this analysis were taken during the 1992-1993 run of the Fermilab 
Tevatron collider and correspond to a total integrated luminosity / Cdt = 6.5±0.4 pb" 1 . The 
dimuon data were collected using a multilevel trigger requiring at least one reconstructed 
muon with transverse momentum p T > 3 GeV/c and at least one reconstructed jet with 
transverse energy Ex > 10 GeV. 

The events are then fully reconstructed offline and subjected to event selection criteria. 
The offline analysis requires two muons with p T > 4 GeV/c and pseudorapidity |?7 M | < 0.8. 
In addition, both muon tracks have to be consistent with originating from the reconstructed 
event vertex and deposit > 1 GeV of energy in the calorimeter. Each muon is also required 



to have an associated jet with E? > 12 GeV within a cone of radius 1Z = y (Ar/) 2 + (A0) 2 < 
0.8. The jet energies are measured using a cone of 1Z — 0.7. Finally, muon candidates in the 
region 80° < M < 110° are excluded due to poor chamber efficiencies near the Main Ring 
beam pipe. 

Further selection criteria are placed on the dimuon candidates to reduce backgrounds to 
bb production. The invariant mass of the dimuons is restricted to the range 6 < m tlfl < 35 
GeV/c 2 . The lower limit removes dimuons resulting from the cascade decay of single b quarks 
and from J /if) resonance decays, while the upper limit reduces the number of dimuons due 
to Z boson decays. An opening space angle requirement of < 165° between the muons is 
also applied to remove contamination from cosmic ray muons. A total of 397 events pass all 
selection criteria. 

The trigger and offline reconstruction efficiencies are determined from Monte Carlo event 
samples. Events generated with the ISAJET [[LI]] Monte Carlo are passed through a GEANT 
|T2| ] simulation of the D0 detector followed by trigger simulation and reconstruction pro- 



grams. Trigger and some offline efficiencies found in this way are crosschecked by using 



appropriate data samples. The overall acceptance times efficiency as a function of the 
higher (leading) muon m the event increases from about 1% at 4 GeV/c to a plateau of 
9% above 15 GeV/c. We define the leading muon in the event as the muon with the greater 
value of pj.. 

In addition to bb production, dimuon events in the invariant mass range of 6-35 GeV/c 2 
can also arise from other sources. These processes include semileptonic decays of cc pairs, 
events in which one or both of the muons are produced by in-flight decays of 7r or K 
mesons, Drell-Yan production, and T resonance decays. Muons from the Drell-Yan process 
and T decays are not expected to have jets associated with them. Monte Carlo estimates 



normalized to the measured Drell-Yan and T cross sections |13j show that less than one event 
is expected to contribute to the final data sample from these two sources. An additional 
source of dimuon events is cosmic ray muons passing through the detector. 

To extract the bb signal, we use a maximum likelihood fit with four different input dis- 
tributions. The input distributions are chosen based on their effectiveness in distinguishing 
between the different sources of dimuon events. We use the transverse momenta of the 
leading and trailing muons relative to their associated jet axes Qo^ 1 ), the fraction of longitu- 
dinal momentum of the jet carried by the leading muon divided by the jet Et (r z ), and the 
reconstructed time of passage (to) of the leading muon track through the muon chambers 
with respect to the beam crossing time. The variable to is used to identify the cosmic ray 
muon background, which is not expected to be in time with the beam crossing. Monte Carlo 
studies show that the variables p!^ 1 and r z help to discriminate between background and bb 
production. For both variables, the jet energy is defined to be the vector sum of the muon 
energy and the jet energy measured in the calorimeter less the expected minimum ionizing 
energy of the muon deposited in the calorimeter. 

The p^ 1 and r z distributions for bb, cc, and b or c plus tt/K decay are modeled using 
the ISAJET Monte Carlo. Each of these samples is processed with a complete detector, 
trigger, and offline simulation. The distributions for b quark decays includes both direct 
(b — > /i) and sequential (b — > c — > fi) decays. The distributions for cc and a c quark plus 
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a 7i or K decay are very similar, so both contributions are fit to the same function. The 
distributions for £ are obtained from two different sources. The t distribution for cosmic 
ray muons is obtained from data collected between collider runs using cosmic ray triggers. 
For beam-produced muons, to is measured using muons from J j if) decays. 

Figure 1 shows the result of the maximum likelihood fit for p^? 1 of the leading muon 
and r z . Included in Fig. 1 are the contributions from each of the major sources of dimuon 
events. The bb contribution to the final data sample is found to be 45.3±5.8%. The other 
fractions fit to the data set consist of b quark plus n/K decay (37.9±5.6%), cc production 
(14.0±3.8%), and cosmic ray muons (2.8±1.6%). From the fit, we obtain the number of bb 
events per bin as a function of pj. of the leading muon and as a function of the difference in 
azimuthal angle between the two muons, A0 M/i . 

The systematic errors on the number of bb events per bin (8%) are estimated by varying 
the input distributions to the maximum likelihood fit within reasonable bounds. As a 
crosscheck of the fitting procedure, we calculate the fraction of events originating from bb 
production using appropriately normalized Monte Carlo samples. Good agreement is found 
between the Monte Carlo calculated fraction and that found from the maximum likelihood 
fit to the data. The fractions agree as a function of both of the leading muon and A0 MM . 



A complete description of the fitting procedure can be found in Ref. [14 



The dimuon cross section originating from bb production is calculated using 

(1) 



dx Ax e(x) J Cdt 

where x is either the pt of the leading muon or A0 MAt . Here, e is the total efficiency, / Cdt is 
the integrated luminosity, N^ 1 is the number of bb events determined from the fit, and f p is 
an unfolding factor to account for smearing caused by the muon momentum resolution. An 
unfolding technique |T3| is used to determine f p . The factor f p varies from 0.78 at low p^. 1 



to 0.93 in the highest pj. 1 bin (/ii is the leading muon in the event) and takes into account 
our invariant mass and pj. requirements. The systematic uncertainty associated with f p is 
found to be p^ dependent and varies from 13% to 22%. 



Figure 2(a) shows the result of the cross section calculation as a function of pt^ 1 for 
4 < Pt < 25 GeV/c, \r]^\ < 0.8, and 6 < < 35 GeV/c 2 . The total systematic error 
is found to be pt^ 1 dependent, ranging from 25% to 31%. This includes uncertainties from 
the trigger efficiency (19%), offline selection efficiency (5%), maximum likelihood fit (8%), 
momentum unfolding (13-22%), and integrated luminosity (5%). 

The theoretical curve of Fig. 2(a) is determined using the HVQJET [[L6[] Monte Carlo. 
HVQJET is an implementation of the NLO calculation of Ref. || (MNR) for bb production. 
It uses the MNR parton level generator and a modified version of ISAJET for hadronization, 
particle decays, and modeling of the underlying event. The particle decays are based on the 
ISAJET implementation |L7| of the CLEO decay tables. In HVQJET the MNR prediction is 
realized by combining parton level events having negative weights with those having positive 
weights and similar topologies. The prediction shown is the NLO calculation and includes 
all four gg, gq, gq, and qq initiated subprocesses with m&(pole mass) = 4.75 GeV/c 2 . The 



MRSR2 |18[ parton distribution functions (PDFs) are used with A 5 = 237 MeV. 

The shaded region in Fig. 2(a) shows the combined systematic and statistical error 
from the HVQJET prediction (1^0%) ■ This error is dominated by the uncertainty associated 
with the MNR prediction and is determined by varying the mass of the b quark between 
4.5 GeV/c 2 and 5.0 GeV/c 2 , and the factorization and renormalization scales, taken to be 
equal, between /x /2 and 2/i , where /i 2 , = m 2 + (p^) 2 . Additional systematic errors include 



those associated with the PDFs (20%), the Peterson fragmentation function [19[] (8%), the B 
meson semileptonic branching fraction (7%), and the muon decay spectrum from B mesons 
(20%). Varying these parameters does not appreciably change the shape of the prediction. 
The Monte Carlo statistical errors are less than 10%. 

To extract the b quark cross section from the dimuon data, we employ a method first 
used by UAl and subsequently used by CDF |]J and D0 0. Since a correlation exists 
between the pr of the muon produced in a b quark decay and the parent b quark pr, cuts 
applied to the muon pt in the data are effectively b quark pt cuts. For a set of kinematic 
cuts, which include cuts on the transverse momentum of the muons, we define p™ m as that 
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value of the b quark px where 90% of the accepted events have b quark transverse momentum 
greater than p™ m . The b quark cross section is then calculated as 

MC 

<p b T>PT n )=^^T)^M^, (2) 

bb— ffj.fi 

where a^{p^) is the measured dimuon cross section of Eq. ([]]) integrated over different 
intervals of p^~ , cr^ 10 is the total Monte Carlo b quark cross section for p T > p™ m (where 
\v b \ < 1.0 and no cut on y b ), and <r^ c is the Monte Carlo cross section for dimuon 

|y 1 a <" bb—>fi,fi, 

production with the same requirements used to select the data set. For each interval of p? , 
p™ m and cr MC are calculated using HVQJET. Combining the uncertainties of the measured 
dimuon cross section with those associated with extracting the b quark cross section, we 
obtain a total systematic uncertainty of 34-38% on the measured b quark cross section. 
The latter uncertainties are associated with b quark fragmentation: Peterson fragmentation 
function; semileptonic branching fraction; and muon decay spectrum with the magnitudes 
noted above. 

Figure 2(b) shows the b quark production cross section for the rapidity range \y b \ < 
1.0 as a function of p™ m . The NLO QCD prediction is computed using Ref. with 
m;,(pole mass) = 4.75 GeV/c 2 and the MRSR2 PDFs. The theoretical uncertainty of Tjs^ 
results from varying the mass of the b quark and the factorization and renormalization scales 
as described above and is dominated by the variation of the scales. The ratio of the data to 
the central NLO QCD prediction is approximately three over the entire p™ in range covered. 

Also shown in Fig. 2(b) is a revised result based on the previous inclusive single muon 
measurement from D0 0. In light of revised B meson decay modes and Monte Carlo 
improvements, the cross section is re-evaluated by using HVQJET to calculate new values 
of cr^ 1 for extraction of the b quark cross section from the measured inclusive single 

muon spectrum. In addition, the high p T inclusive muon data (pj. >12 GeV/c) are excluded 
due to large uncertainties in the cosmic ray muon background subtractions. The resulting 
increase in the b quark cross section is primarily caused by the new B meson decay modes 
and lower semileptonic branching fractions ||17|| . The re-evaluated cross section supersedes 
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that of Ref. @. The tabulated data for the dimuon and inclusive single muon data sets can 
be found in Tables I and II. 

The differential bb cross section, da^ / dA0 MM , gives further information on the underlying 
QCD production mechanisms. The azimuthal opening angle between b and b quarks (or 
between their decay muons) is sensitive to the contributing production mechanisms. These 
contributions are the leading order (LO) subprocess, flavor creation, and the next-to-leading 
order subprocesses, gluon splitting and flavor excitation. There are also contributions from 
interference terms. 

The cross section da^ - / dA^^ is shown in Fig. 3. Also shown are the LO and NLO QCD 
predictions which are determined using HVQJET and include all subprocesses. The grey band 
around the NLO prediction shows the combined statistical and systematic errors associated 
with the prediction, which is ^50% as detailed above. The data again show an excess above 
the NLO QCD prediction but agree with the overall shape. The agreement in shape is 
consistent with the presence of NLO subprocesses since the LO prediction, which contains 
the smearing from the b — > B — > /x fragmentation and decay chain, does not describe the 
data. 

In conclusion, we have measured the b quark production cross section and the bb az- 
imuthal angle correlations using dimuons to tag the presence of b quarks. These measure- 
ments, as well as the revised inclusive single muon measurement, are found to agree in shape 
with the NLO QCD calculation of heavy flavor production but lie above the central values 
of these predictions. 
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FIGURE CAPTIONS 



Figure 1 The results of the maximum likelihood fit to the data for (a) p™ 1 of the leading 
muon and (b) r z . Also included are the curves showing the contribution from each process 
to the dimuon sample. 

Figure 2 (a) The unfolded leading muon p T spectrum for bb production compared to the 
predicted spectrum (see text) where the data errors are statistical (inner) and total (outer) 
and the Monte Carlo errors are total (shaded band); (b) the b quark production cross section 
for \y b \ < 1.0 compared with the revised inclusive single muon results and the NLO QCD 
prediction. The error bars on the data represent the total error. The theoretical uncertainty 
shows the uncertainty associated with the factorization and renormalization scales and the 
b quark mass. Also shown are the inclusive single muon data from CDF [|TJ. 

Figure 3 The A0 W spectrum for bb production compared to the predicted spectrum (see 
text). The errors on the data are statistical and total. The solid histogram shows the 
NLO prediction with the grey band indicating the total uncertainty. Also shown is the LO 
prediction (dotted histogram) with the statistical error only. 
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Figure 3 



bb^ /.w, X 

□ Data 

— HVQJET 
HVQJET LO 



[] 



[] 



[] 




20 40 60 80 100 120 140 160 18 

[Degrees] 



19 



TABLE CAPTIONS 



Table 1 Cross sections for bb — > fifi production. 



Table 2 Results for the b quark production cross section for \y b \ < 1.0. 
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TABLES 



Table 1 



Ml 

Pt 


daw/dpr 


Stat Error 


Syst Error 


[GeV/c] 


[nb/GeV/c] 


[nb/GeV/c] 


[nb/GeV/c] 


4-5 


0.30 


0.095 


0.092 


5-7 


0.21 


0.018 


0.054 


7-10 


0.064 


0.0037 


0.016 


10-15 


0.018 


0.0012 


0.0046 


15-25 


0.0025 


0.00034 


0.00067 
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Table 2 



Pt 


a b 


Total Error 


Data Source 


[GeV/c] 


[fib] 


[Mb] 




6 


8.76 


2.80 


Single Muon 


7 


7.78 


1.95 


Single Muon 


8 


5.96 


2.29 


Dimuon 


9 


4.85 


1.22 


Single Muon 


11 


3.10 


1.08 


Dimuon 


12 


2.28 


0.58 


Single Muon 


13 


1.71 


0.59 


Dimuon 


15 


1.16 


0.30 


Single Muon 


18 


0.61 


0.21 


Dimuon 


25 


0.17 


0.06 


Dimuon 
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